Near-threshold Photoproduction of (f) Mesons from Deuterium 
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We report the first, kinematically-complete measurement of the differential cross section of <fi- 
meson photoproduction from deuterium near the production threshold for a proton using the CLAS 
detector and a tagged-photon beam in Hall B at Jefferson Lab. The measurement was carried out 
by a triple coincidence detection of a proton, K + and K~ near the theoretical production threshold 
of 1.57 GeV. The extracted differential cross sections 4z for the initial photon energy range of 1.65- 
1.75 GeV are consistent with predictions based on a quasifree mechanism. Our finding is different 
from recent LEPS results on </>-meson photoproduction from deuterium in a similar incident photon 
energy range, but in a different momentum transfer region. 

PACS numbers: 13.60.Le, 24.85. +p, 25.10.+S, 25.20.-x 



The subject of medium modification of properties of 
(light) vector mesons such as their masses and widths 
has been an important subject in nuclear physics, par- 
ticularly in the context of partial restoration of chiral 
symmetry [l[. Among light vector mesons, the <fi me- 
son is considered an excellent candidate for the study of 
medium modification effects due to its narrow, vacuum 
width of 4.3 MeV/c 2 . The first evidence for in-medium 
modification of the (fi meson at normal nuclear density 
from hadronic reactions was reported by the KEK-PS 
E325 Collaboration 0], in which a mass shift of 3.4% and 
a width increase of a factor of 3.6 for (fi momenta around 
1 GcV/c from reaction p + A —> p,u!,<fi + X — >• e + e~ + X' 
were deduced. 

Studies of </>-meson production from nuclei have gener- 
ated new interest in the last a few years motivated largely 
by a number of new and interesting results reported by 
different experiments. A number of experiments reported 
larger values of (fi-N total cross section extracted from 
photoproduction of 0-meson from nuclear targets than 
an upper limit of n ~ 11 mb Q obtained using the 
0-meson photoproduction data on the proton and the 
vector meson dominance (VMD) model |4|. This upper 
limit is in agreement with the estimate from the addi- 
tive quark model The LEPS collaboration re- 
ported the A-dependence of </)-meson photoproduction 
from Li, C, Al, and Cu nuclear targets, and a larger 
value of (j07v — 35±^ mb was obtained, correspond- 
ing to a collisionally broadened width of 35 MeV, much 
larger than that reported by the KEK-PS E325 Collabo- 
ration The CLAS Collaboration also reported larger 
values of u^n from a combined analysis [7[ of the co- 
herent Q and incoherent 0-meson photoproduction from 
deuterium and a value larger than 20 mb is favored for 
t7(f>N- Another recent CLAS experiment [9| reported a 
measured value in the range of 16-70 mb for u^n from 
■yA -> <fiX -> e + e~X. More recently, the Anke Col- 
laboration at COSY reported [13] a nuclear dependence 



of ^4°- 56 ±o.03 f rom pA —j. (fix in the <fi meson momen- 
tum region of 0.6 - 1.6 GcV/c. The comparison of the 
data with model calculations suggests an in-medium <fi 
width about an order of magnitude larger than its vac- 
uum value, which is consistent with a larger value of 
o~4>n @] taking into account the experimental uncertain- 
ties of both experiments. 

In order to determine the medium modification of the 
properties of <fi unambiguously, understanding reaction 
mechanisms of </>-meson production from nuclear targets 
is crucial because nuclear targets are necessary in ex- 
periments studying the medium effect. The deuteron, a 
loosely bound and the simplest nuclear system, provides 
an excellent testing ground to study the reaction mech- 
anisms. Such a study has been reported [IH recently 
by the LEPS collaboration on the measurement of the 
incoherent "fd — > cfipn photoproduction near threshold 
at forward angles. The nuclear transparency ratio for 
(fi has been extracted, which shows a large suppression 
consistent with the nuclear dependence reported Q by 
the same collaboration previously using 0-meson photo- 
production data from heavier nuclear targets. However, 
a recent theoretical study 12j of 0-mcson photoproduc- 
tion from a deuterium target shows that the contribution 
from rescattering at the kinematics of the LEPS experi- 
ment is small compared with the cross section of jp — > cfip 
in free space. 

Another important aspect of </>-meson production from 
nuclear targets near threshold is related to experimen- 
tal investigations of the existence of a <f> — N bound 
state. The QCD van der Waals interaction, mediated by 
multi-gluon exchanges, is dominant when two interacting 
color-singlet hadrons have no common valence quarks. 
Bound states of i] c with 3 He and heavier nuclei were pre- 
dicted [HQ. Similarly, one expects that the attractive 
QCD van der Waals force also dominates the (fi-N inter- 
action. A bound state of (fi-N was found to be possible by 
Gao, Lee, and Marinov [15j|. Such a bound state was also 
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predicted by Huang, Zhang, and Yu using a chiral SU(3) 
quark model and the extended chiral SU(3) quark model 
by solving the Resonant Group Method (RGM) equa- 
tion fl6| . Subthreshold production kinematics have been 
proposed 
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as advantageous for the search of <fi — N 
bound states from nuclear targets. 

In this paper, we report on the first kincmatically- 
complete measurement of differential cross sections from 
near-threshold production of 4> mesons from a deuterium 
target by a triple coincidence detection of K + , K~ , and 
proton. The measurement was carried out using the 
CLAS detector [l?} at Jefferson Lab. The incident pho- 
ton energy range used in this analysis is 1.65 - 1.75 GeV, 
which is above the (/>-meson photoproduction threshold 
^fithres. _ 2 57 GeV) from a free proton target. How- 
ever, due to the requirement of a triple-coincidence de- 
tection and the imperfect acceptance of the detector at 
forward angles, the reconstructed <f> event in this anal- 
ysis originated mostly from photoproduction on a high- 
momentum proton inside the deuteron, and is below the 
CLAS production threshold for -fp — >• (f>p — >• K + K~p. 
Such production therefore images subthreshold produc- 
tion of (/>-meson from nuclear targets and as such is im- 
portant for future experimental search for a </>-N bound 
state. 

High statistics data were collected during the CLAS 
glO running period [l8j from a 24-cm-long liquid- 
deuterium target. A tagged-photon beam was used, 
which was generated by a 3.8-GeV electron beam inci- 
dent on a gold radiator with a thickness of 10~ 4 radi- 
ation lengths. The photon flux was measured by the 
Hall B photon-tagging system [lj|. Two settings of the 
CLAS magnetic field were used during the experiment, 
corresponding to a low-field setting (with a toroidal mag- 
net current 1=2250 A) for better forward-angle coverage, 
and a high-field setting (1=3375 A) for better momen- 
tum resolution. The reaction d("/, <f>p)n was measured by 
detecting kaons from the (/)-meson decay (4> — > K + K~, 
branching ratio about 0.5), using the same data set as in 
Refs. 0,1,0]. 

The K + , K~, and the proton were selected based on 
the particle charge, momentum, and time-of-flight in- 
formation. The reaction di^y, 4>p)n was identified in the 
missing mass squared distribution by applying a ±3cr cut 
on the missing neutron peak. The energy threshold for 
the 7A — > (f)N reaction is 1.57 GeV. However, due to the 
minimum detection threshold for charged particles, the 
CLAS acceptance determined threshold is around 1.75 
GeV. This is demonstrated in our analysis as no 4> events 
at incident photon energies below 1.75 GeV can be iden- 
tified from the glO hydrogen data set, which was taken 
during the glO running period for calibration purposes. 
This finding is consistent with the hydrogen results from 
the CLAS gll |2l[ high statistics data set. Table Q] sum- 
marizes the differences between the glO and gll exper- 
imental settings. Fig. [1] (left panel) shows the invariant 



Run 


Magnet 


Target Target Target 


Accumulated 


Period 


Current 


Material Length Position 


Flux 


gio 


2250 (A) 


LD 2 24 cm -25 cm 


1.266X10 12 


glO 


2250 (A) 


LH 2 24 cm -25 cm 


7.508xl0 10 


gll 


1930 (A) 


LH 2 40 cm -10 cm 


4.316xl0 12 



TABLE I: A comparison of the glO and gll experimental 
settings. The accumulated photon flux is for the Ej range of 
1.65-1.75 GeV. The "-10 cm" means that the target center of 
the gll hydrogen target is 10 cm upstream from the nominal 
center of CLAS. 



mass distribution of the K + K~ before the acceptance 
correction from the glO deuterium data set, where the <fi 
peak is clearly visible. Also shown in Fig. [T] (right panel) 
is the corresponding spectrum from the gll hydrogen 
data set normalized to the glO integrated luminosity for 
nucleons. The yield of the gll hydrogen is ~ 1.5 events 
per 2.5 MeV around the reconstructed cj> meson mass, 
which is strongly suppressed compared to that of the glO 
deuteron (~ 23 events per 2.5 MeV) due to the energy 
threshold of producing a (f> meson on a nucleon at the 
kinematic settings of glO and gll. Therefore, the pho- 
ton energy range used to extract the near-threshold cross 
section for 0-meson photoproduction from deuterium is 
between 1.65 GeV and 1.75 GeV in this work. The chosen 
photon energy range was further confirmed by simulating 
the 7 + p — >• p + 4> process for the glO configuration. 

Once the reaction d( r y,pK + K~)n was identified, the 
number of <f> mesons was obtained by subtracting the 
background under the <f> peak (invariant mass spectrum 
of the K + and K~) in the ±3cr region (see Fig. [T]). The 
K + K~ invariant mass distribution was fitted using a 
Breit-Wigner function convoluted with the experimen- 
tal resolution, plus a function to model the background 
in each kinematic bin. The experimental resolution on 
the missing mass ranged from 1.2 to 1.7 MeV for differ- 
ent t bins. They were obtained by fitting the invariant 
mass distribution of the Monte-Carlo simulation. The 
background shape was assumed to be [1] : 

f(x) = a^x 2 - (2M K ) 2 + b (x 2 - (2M K ) 2 ) for x > 2M K 
f(x) = for x < 2M K , (1) 

where x is the invariant mass of the K + K~, Mk is the 
kaon mass, and a and b are fitting parameters. Such a fit 
was performed separately for each t bin, and the t depen- 
dence of the background was effectively included in the 
fitting parameters a and b. In addition, the background 
was also fitted to a straight line. The results from fitting 
these two shapes were compared in order to estimate the 
systematic uncertainties due to the subtraction of the 
background. 

GEANT3 Monte-Carlo (MC) simulations were carried 
out to model detector efficiencies and resolutions for this 
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FIG. 1: The K K invariant mass distribution from the 
CLAS glO deuterium data for £ 7 = 1.65-1.75 GeV (left 
panel). The a;-bin size is 2.5 MeV. The corresponding dis- 
tribution from the CLAS gll hydrogen data set in the same 
photon energy range is shown in the right panel. 
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reaction channel. A quasifree event generator was used 
for the near-threshold kinematics. It generated pK + K~ 
three-body events with a random photon energy based on 
the photon energy distribution of the data in the energy 
range of interest. The initial momentum of the nucleon 
inside the deuteron was chosen using the Bonn potential 
wavefunction |22| . The spectator nucleon was assumed to 
be on-shell, whereas the struck nucleon was assumed to 
be off-shell before absorbing a photon. Isospin symmetry 
was assumed for the (/)-meson photoproduction from the 
nucleon. The events were then checked to ensure energy 
conservation. The MC events were generated based on 
a Brcit-Wigncr shape of the resonance centered at the <\> 
mass of 1.019 GeV 2 with a full width at half maximum 
(FWHM) of T=4.26 MeV. The 4> meson decay angular 
distribution and cross section are based on the gll hy- 
drogen data Q ■ 

The ^>-meson differential photoproduction cross section 
on a hydrogen target (% vs. t — to) wa s obtained using a 
fit to the gll data in £L = 1.625 — 3.775 GeV range. Here 
t is the four-momentum transfer squared, — -P 7 ) 2 , 
and to is the maximum t value for a given photon energy 
because t is negative. In addition, the event generator 
included the TV — TV and — TV final-state interactions 
(FSIs). The Jost function approach |24| was used for 
the TV - TV FSI. The </> - TV FSI, which was assumed to 
be incoherent from the original 0-mcson photoproduc- 
tion process, was modeled based on the vector meson 
dominance (VMD) model, in which the t-dependence of 
the 4> — TV elastic scattering cross section is the same as 
that of </>-meson photoproduction. A fitting procedure, 
in which the strength of TV — TV and cf> — TV FSIs were 
obtained, was then used to optimize the 0-meson pho- 
toproduction model so that the resulting Monte-Carlo 
distributions match those of the data. Fig. [2] shows the 
comparison in the missing momentum distribution for 
the data (solid triangles) and the MC (solid circles). 

MC-generated events were used as input to the 



FIG. 2: Missing momentum distributions are shown for data 
(solid triangles) and MC (solid circles), where the MC results 
are normalized to the luminosity of the data. Both (/> — TV and 
TV - TV FSI are included in the Monte Carlo. 



GEANT3-based CLAS simulation |25|. They were then 
reconstructed using the same algorithm as was used for 
the data. The acceptance was obtained by the ratio of 
the number of events that passed the analysis cuts to the 
number of generated <fr events. The average differential 
cross sections were extracted by dividing the normalized 
yield by the acceptance. The differential cross sections 
were then bin-centered at fixed t values with a finite bin- 
ning correction. 

Several sources contribute to the overall systematic un- 
certainty in the differential cross section. The systematic 
uncertainties associated with particle identification and 
the missing mass cut were 4.2 - 12.9% and 1.4 - 10.9%, 
respectively. These values were determined by varying 
the corresponding cuts by ±10% in each t bin. The an- 
gular distributions of the </>-meson's decay products in 
its rest frame and the cos(# c . m .) distribution of the <p 
meson were uncertain to within 10% and 5% 0, l26j |. 
respectively, leading to 5.2-13.2% systematic uncertain- 
ties. The background obtained from the non-linear back- 
ground shape was on average 8% smaller than that from 
the linear background. A conservative 8% systematic un- 
certainty is assigned for the background subtraction pro- 
cedure. The systematic uncertainties due to the effect of 
FSIs were obtained by varying the fitted strengths of the 
TV — TV FSI and the - TV FSI by 30% and 50%, respec- 
tively. The systematic errors vary from 4% to 17% for 
different t bins. The uncertainty in the photon flux was 
20l 27 1 . The uncertainty of the bin-centering correc- 



5% 

tion was assumed to be 30% of the size of the correction 
based on knowledge of the CLAS acceptance and the in- 
put cross section model, which is obtained from the gll 
data [23[ . The absolute size of bin-centering corrections 
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FIG. 3: (Color online) <j> photoproduction differential cross 
sections from deuterium are plotted as a function of \t — to\. 
The inner error bars are the statistical uncertainties, and the 
outer error bars are the quadrature sum of systematic and 
statistical uncertainties. 
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stat. 


sys. 


(GeV) 


(GeV 2 ) (/ib/GeV 2 ) 


uncer. 


uncer. 


1.65-1.75 


0.509 0.31 


0.084 


0.094 


1.65-1.75 


0.887 0.20 


0.049 


0.037 


1.65-1.75 


0.924 0.13 


0.066 


0.041 



TABLE II: Tabulated results on the <^>-meson photoproduc- 
tion from deuterium for an _E 7 range of 1.65 to 1.75 GeV. 



was between 1% and 10%. Combined in quadrature, the 
overall systematic errors vary from 18% — 32% depending 
on the kinematics. 

The <^-meson photoproduction cross sections for the 
deuteron are tabulated in Table HI] and are plotted as a 
function of \t — to\ in Fig. [3] for a photon energy range of 
1.65-1.75 GeV. The solid circles arc results obtained from 
the CLAS glO low-field setting, whereas the solid square 
is the result from the high-field setting. The quasifree 
calculation is also plotted for comparison together with 
its uncertainty (shown as a band in Fig. [3J . This sim- 
ple calculation is based on a quasifree picture with the 
0-meson differential photoproduction cross section from 
the proton which is based on the gll data [23J. The prin- 
ciple of this calculation is the same as that for the event 
generator used in the MC. The systematic uncertainty for 
this calculation is about 10% due to the uncertainty in 
the input cross section. The extracted differential cross 
section for 0-meson photoproduction is consistent with 
the quasifree calculation within uncertainties. Our find- 
ing is consistent with theoretical calculations presented 
in Ref. [l2| . However, they are not consistent with the 



recent LEPS results ll( from deuterim in a similar pho- 
ton energy range, but in a smaller momentum transfer 
region. The corresponding \t — *o| range for the LEPS 
data in a photon energy range of 1.65 to 1.75 GeV is less 
than 0.2 GeV 2 . Future studies both in experiment and 
in theory are important to clarify the situation. 

In summary, we have extracted for the first time the 
differential cross section on 0-meson photoproduction 
from deuterium from a kinematically-completc measure- 
ment below the production threshold for the proton ac- 
cessible on CLAS. The chosen incident photon energy 
range is 1.65-1.75 GeV, which is near the 1.57 GeV pro- 
duction threshold for protons. Our extracted differential 
cross sections are in agreement with predictions from a 
simple quasifree picture. This finding is consistent with 
the recent theoretical study of this reaction [l2j , though 
inconsistent with the recent LEPS results in a somewhat 
different kinematic region. Further, our data provide 
information on 0-meson subthreshold photoproduction 
cross section from a deuterium target. Although heavier 
nuclear targets will be ideal for future dedicated searches 
for a 0-N bound state, the extracted cross sections from 
deuterium reported in this paper will help provide infor- 
mation on the expected production rate of the </>-N bound 
state from heavier nuclear targets. 
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